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COLLISION INTEGRALS FOR A MODIFIED STOCKMAYER POTENTIAL 

By Eugene C. Itean, Alan R. Glueck, 
and Roger A* Svehla 


SUMMARY 

Collision integrals we re calculated for the modified Stockmayer 

potential E(r) = 4e[(cr/r)-^ - (a/r)^ - D(o/r)^], which maybe applied 
to polar molecules. It was assumed that the colliding molecules main- 
tain their same relative orientation during the encounter. Calculations 
of the integrals were made for a large reduced temperature range and for 
a range of 5 from 0 to 10. The results agree with other work on non- 
polar interactions (6 = 0). However, for polar interactions, the only 
previously published calculations have been found to be in error and do 
not agree with this work. 

Assuming that the molecules interact as alined dipoles of maximum 
attraction, values for a, e, and & were determined for various polar 
molecules by a least squares fit of experimental viscosity data. Satis- 
factory results were obtained for slightly polar molecules, but not for 
more highly polar molecules such as NH 3 or H 2 0. Therefore, it appears 

that the assumed model of molecules interacting at all times as alined 
dipoles of maximum attraction is not satisfactory for estimating trans- 
port properties of polar molecules. 


INTRODUCTION 

Coefficients of viscosity, thermal conductivity, and diffusion are 
needed in heat- and mass -transfer calculations. Equations for calcula- 
tion of these transport properties have been developed from kinetic 
theory in terms of collision integrals, quantities that describe the 
interaction between colliding molecules. When these integrals axe known 
it is possible to predict transport properties at elevated temperatures. 
Values of these integrals (ref. 1 , pp. 1126-1180) have been calculated 
assuming various interaction potentials. However, these integrals are 
specifically for nonpolar molecules. In many cases a polar gas such as 
H 2 0 , NH 3 , or HC1 is of interest, and the predictions would be in error 



? 


if these integrals were used, "because the pol.ir character of the gas is 
ignored. 


The Lennard-Jones potential of interaction for spherically symmetric 
nonpolar molecules (ref. 1, p. 32), given in equation (l), is a well- 
known potential that has been successfully used for correlating transport 
properties of many nonpolar gases: 


E(r) 



where E(r) is the potential energy of interaction, r is the inter- 
molecular distance between colliding molecules, € the maximum energy of 
attraction, and a the value of r where the potential energy of in- 
teraction is zero. For low-velocity encounte rs, a can be considered 
the collision diameter of the molecule. (All symbols are defined in 
appendix A. ) These relations are shown in figure 1. 


For polar molecules, an equation similar to equation (l) has been 
proposed by Stockmayer (ref. 2), given in a medified form in equations 
(2) and (3): 


E(r) = 4e 


ir - $ 




;( 0 2 j?) 


g(0q, 02^^) = 2 cos 0 q_ cos 0 2 - sin sin cos 


( 2 ) 

( 3 ) 


Stockmayer actually proposed this equation using an arbitrary exponent 
S, instead of 12 as given in equation (2). Bat written in the preceding 
form it can be considered a Lennard-Jones potential modified to include 
the forces between two point dipoles. The angles 0q and 0£ are the 

angles of inclination of the dipoles with the intermole cular axis, p is 
the dipole moment of the molecule, and 9 is the azimuthal angle between 
the dipoles. This is shown in figure 2. However, the molecular con- 
stants a and e do not have quite the same significance as in equa- 
tion (l). They now represent constants that would be obtained from the 
interaction of a polar and a nonpolar molecule. 

Define the parameter & as follows: 

& = — ti — — g( 0^_j 
4 € 0 ° 

Using this definition, equation (2) may be rewritten as follows: 

i?f - (d 6 - < 


E(r) = 4e 


(5) 
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Equation (5) is the form of the equation that Krieger (ref. 3) used to 
obtain collision integrals for polar molecules. In order to simplify 
calculations, Krieger suggested letting g(0q,0 2 , cp) in equation (4) equal 

+2, corresponding to alined dipoles of maximum attraction. This assumes 
the dipoles have sufficient time to orient before collision. For this 
assumption, equation (4) becomes 



For the purpose of calculating the collision Integrals using equation 
(5), no specific orientation of the colliding molecules need be assumed. 
If it is only assumed that the molecules maintain their same relative 


orientation during the encounter, then g (0-^ 0 2 ,cp) becomes a constant, 


and it is possible to assign constant values of 5. Krieger performed 
his calculations for positive 6 values. Since a wider range of param- 
eters was desired, and because of a lack of smoothness of Krieger' s re- 
sults, collision integrals for positive values of & were recalculated. 


The calculation of the integrals requires three integrations. The 
first Is to obtain the angles of deflection, the second the collision 
cross sections, and the third the collision integrals. These three in- 
tegrals (ref. 1, pp. 525 - 527), based on equations given by Chapman and 
Cowling (ref. 4), are given in the next section. 


The collision integrals evaluated were the and the at 1 .? 1 )* 

integrals. These are the ones used In evaluating first approximations 
to the coefficients of viscosity, thermal conductivity, and diffusion of 
pure gases. The significance of the superscripts is indicated in the 

following section. Other integrals such as fit 1 - 2 )*, q(X, 3)* or a (2,3)* 
could be calculated in a similar manner. These are used In higher approx- 
imations or for mixtures. But the approximate nature of the potential 
assumed does not warrant their use. 


The transport property equations for pure gases (ref. 1, pp. 528- 
539) are given here for convenience: 


huxio 7 = 

[DjxiO 7 = 0 i P02628QVt5/M 



(7) 

( 8 ) 
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where 

D self -diffusion coefficient, cm 2 / sec 

k Boltzmann's constant, ergs/°I< 

M gas molecular weight, g/mole 

p pr e s s ur e y atm 

T temperature ^ °K 

■q viscosity, g/(cm)(sec) 

o collision diameter 9 angstroms 

^(l,l)*^(2,2)* collision integrals from table II evaluated at reduced 

temper at i ore, T* = kT/e 

No equation is given for thermal conductivity, because, for polyatomic 
molecules, the equation is complicated by the interconversion of trans- 
lational and internal energy. In polar molecules, additional con lg- 
urational effects are involved (refs. 5 and 6). 


DERIVATION AND EVALUATION OF COLLISION INTEGRALS 
Equations of Mot i or 

The general equations of motion for a tvo-body system in polar co- 
ordinates, consisting of two identical colliding molecules A and B, and 
describing the motion of A relative to B, are as follows: 


- m(r 2 0 2 + r 2 ) + E(r) = j mg 2 


( 10 ) 


where e is the relative velocity before the encounter at r - °°, 
b is the distance of nearest approach in the absence of any interactions 
m is the mass of each particle, r is the i itermoleculax separation, 
and E(r) is the potential function given in equation (5). These re 
lations are Illustrated in figure 3. 


If time is eliminated between equations (9) and (10), 


equation is 




-1 



the resulting 


(ii) 
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Inf cgro/tirig from r = a “to r = wh6rs 9- is "the cListojics of closes! 

approach, results in the expression for the angle Q^, the angle of 6 
for r = 00 : _ _ 

^2 4r^E(r)l"^' ^ 




00 P 

.*.»> - / 


— - 1 


mg 2 b 2 


1 ' 


( 12 ) 


1 

■q 


Since dr/d.6 = 0 at r = a, equation (ll) becomes 

a£ _ ± 4a 2 E(a) = Q 

b 2 mg 2 b 2 


(13) 


which determines the lower limit a. Knowing 0 m as a function of g 
and b, the cross section for transport can be found from the relation 
(ref. 1, p. 525) 


Q (l) (g,6) = 2* f (1 - cos^) X)b 

* / 0 


db 


(14) 


(15) 


where l is a positive integer, and 

X = * " 2^ 

as shown in figure 3. 

In this work the only cross sections of interest are Q(l)(g) and 
Q^ 2 )(g). Equation (14) then becomes 


Q vl \gj5) = 2it j 

f (l + cos 20,0b db 

0 

(16a) 

Q^(g,&) » 2 it 

f%CO 

j (1 - cos 2 2Q m )b db 

(16b) 


From the cross sections, the final collision integrals can be computed 
(ref. 1, p. 525) by 


n(l> s )(T, 6 ) = y* e~ r2 r 2s+3 Q^^(g) IT 


(17) 


where s is a positive integer, and 


r 2 = mi 

1 4kT 


The collision integrals calculated were ft(l>l)(T) and o( 2 ^ 2 )(T). 


( 18 ) 
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Introduction of Dimensionless Parameters 

As an aid to computation, the variables can be put in dimensionless 
form, utilizing the characteristic quantities e and o Define the 
following dimensionless quantities: 


r* ^ r/o 
b* = b /a 


E* = E(r*)/e = 4(r 


*-12 _ r *-6 _ 5 r *- 3 ^ 


Letting 


g*2 _ rn g2^4 € (2; 

T* = ET/€ (2: 

and using the values of Cp' and for a rigid sphere (ref. 1 

p. 525) 


qO ) 

rigid sphere 


r . 1 1 + ("pH 

L 2 1+1 j 


fl (M) 

rigid sphere 


= t/H 

▼ Jtra 


t s + l)! 

2 


'rigid sphere 


then equations (12), (14), and (17) become in reduced form (ref. 1, 
p. 527) 


Sm(g*>* 6) 


= / dr* /r 

*4* r * V 3 


]-* 2 e* \ ly ^ 2 


g* 2 b* 2 


Q ( ^*(s*,6) = 


" I 

rigid sphere 1 1 - 


1 1 ± (-D 1 

2 l + l 


[l - COS U)x)b< 


n (l,s)*( T * g) = 


rigid sphere 


(s + l) !T* S+ ^ 


r e -(g* 2 /T*) „(l)* (g . )dg . 

Jr\ 
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In particular. 


q (!)' 


/OO 

(l + cos 20 ra )b* db* 


(29a) 


Q (2) *(gl&) =3 f (1 - cos2 2 ^)b* db* 

J o 

/ OO 

e-(«* Z /T*) g#5 <j(U V)d«* 


(29b) 

(30a) 


0<M>V,#> = r^j f e'^/TP g ,7 Q (0 V)dg* 

5T y n 


(30b) 


Integration Technique 

Angle of deflection integral * - The angle of deflection integral 


^(g*,t>*,5) 



(26) 


In order to simplify the numerical integration, the change of variable 
p = l/r is used in equation (26) to give 


0m(g*,b*,S) 



*2 


11/2 


b*2p2 + 


,*2 


(-p 12 + p 6 + 6p 5 ) 


dp (31) 


where A(= l/a*) is the smallest positive root of the denominator of 
the integrand of equation (3l): 


~~2 ( “P 12 + P 6 + 5p3) - b + 1 = 0 (32) 

§ 

Equation (32) has either three or one positive roots (by Descartes' 
rule of signs), so that finding the smallest root proves to be a problem 
in some cases. Figure 4 shows a typical example of this function for 
g =0.5, 5=1^ and various values of b*. 
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To find the proper root of equation (32), an iterative procedure ^ 
(Uspensky* s method, ref. 7) is used with the initial estimate A — l/h * 
This is a good estimate for sufficiently large h , and also for all 
■ft* < since an examination of equation (32) shows that A is ^greater 
than unity for all b* 1. However , in certain cases where h i° 
not much greater than 1, a minimum above the p-axis exists for equa- 
tion (32) in the vicinity of the initial estimate p = l/b (e.g., 

-p* - 3 in. fig- 4). In this situation Uspensky’s method does not con- 
verge but oscillates about the minimum. This situation can be detected 
by the fact that the first derivative becomes positive during this 
oscillation. When this occurs a new estimate of A > 1 is. made so that 
the first derivative of equation (32) is negative, and the iterative 
procedure is continued. This root-finding procedure has been successful 
In all cases attempted in this work. 

Once the proper root has been founds the angle of deflection inte- 
gral can be evaluated. It should be noted that a pole exists at the 
upper limit A. However, it is a half-order pole if A is a single root, 
and can be handled by fitting the function to a polynomial of the form 

(aQ + a x x + . . • + a n x n )/x 1 / 2 = f(x) in the -icinity of the pole. The 

remainder of the integral is well "behaved and is accurately evaluated 
using the Gaussian numerical integration procedure (ref. 8). 

ppf’pvTlnv to figure 1. it is seen that a double root of equation 
(32) if possible (in this Sample at b* - bg, where b£ s 3.6). The 

integral now no longer has a half -order pole at A, but rather a pole 
of order 1, which leads to an integral whose value approaches ». Thus, 
for values of b* near b£, the molecules orbit around each other an 

indefinite number of times before separating. The occurrence of a 
double root is not possible for all values of the parameter g . For 
every value of & there exists a value Sq(^)j such that for all 

0 < g* ~ g* there exists a b^(g*, &), where aquation (32) has a double 
root, and the phenomenon of "orbiting” occurs. For all g > gQ no 
positive value of b* exists such that equation (32) has a double root, 
and orbiting cannot occur. Values of g£ 2 (&) are given in table I. 

Cross-section integrals. - The cross-section integrals, equations 
(29a) and (29b), can be easily evaluated for g* > gQ(b) by dividing 
the integral into two parts: 

(a) 0 to b^ 

(b) b* to co 
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where for all b* > bjj, it has been found empirically that 6^ can be 
well approximated by 


That is to say, 0 m approaches ji/ 2 asymptotically as b* 


(33) 


The first region is evaluated numerically using the Gaussian inte- 
gration procedure. The second region is evaluated by substituting equa- 
tion (33) for 0 m , expanding cos 2 0 m in a Maclaurin series, and inte- 
grating equations (29) analytically. Dropping all terms after the first 
nonvanishing term gives: 


/ 


R 


/: 


(l + cos 26^* db* s i 


(1 - cos 2 20 m )b* db* s 



46 


e* 2 b# 2 
i S d r 


(34a) 


(34b) 


R 


When g* < g^(6), orbiting occurs and the curve ^ against b* has a 

singularity at 1 q. For this situation the integral was broken into 
five parts as in reference 9 (see fig. 5): 


(a) 

0 

to 

'’M 

Cb) 


to 

b 0* 

(c) 

*0 

to 

b* 

(a) 


to 


(e) 

bR 

to 

00 

Regions 

(a, 

) and (d) are evaluated numerically as was done in the 


first case. Also, region (e) is evaluated as before by using equations 
(34). The regions (b) and (c), which are in the neighborhood of the 
singularity at b^, are evaluated by curve-fitting 6^ against b* 2 

by an empirical equation of the form 

0ra = a ° + , *2 1 . *2 (35) 
b Q - b 


where aQ and ap are constants. 



10 


If this substitution is made for 0^ in equations (29), the fol- 
lowing results are obtained by integrating analytically: 



* r - r / 
0 . _ /v* 2 _ -k*2^ | ■ 


(1 + cos 28)d(b*^) = (bJ J - b M ) + cos 2% 

- 2 (% - a 0 ) |(cos 2a 0 )[si(2% - 2ao) - + (sin 2ao)Ci(2% - 2ao)|^ 


(36a) 



- 4 (% ■ a o) 


? (bj 2 - b^ 2 ) / 

29)d(b* 2 ) = g ( 1 - cos 4% 

|(cos 4a o )Jsi(40 M - 4a Q ) - |j + (sin 4a o )Ci(40 M 



(36b) 


where 


and 


Si(x) = 


Ci(x) = 



Similar results are obtained for the region :rom bp to 



Collision integrals . - The collision integrals are given by equa- 
tions (30a) and (30b). The final integration that obtains the collision 
integral is divided into two parts i 

(a) 0 to gQ 

(*>) gj to g^ 

The integral over both parts is evaluated numerically using 
Gaussian integration, the only difference being in the manner in which 
the cross sections q( 0* are calculated. This has already been dis- 
cussed. The integration is terminated at seme g* = g* where the in- 
from g^ to °° is negligible compared with the total integral. 

For all T* S. 512, g^ is less than 120. 


I 


TfiA 



tact 
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DISCUSSION OF RESULTS 
ow UltS ° f the calcul ation of the collision integrals fl(l,l)' 


are given in table II. The values extend over 


a larKe 


T - 512, and ten values 
selected for ease in inter- 


and 

reduced temperature range from T* = 0. 25 to 
Of 5 from 0 to 10. The T* intervals were 
polation and for comparison with other work. 

Results of this paper, Hirs chf elder ' s values (ref. l), Krieger's 

results (ref. 3), and Rowlinson's values (ref. 10) of ft(2,2)* for 

6 = 0 showed agreement. Moreover, Hirs chf elder’s qUU)* fnT . * _ n 

agreed with this work. However, for 6 > 0 the results of Krieger 2nd 

KrllZ?s vaC* a « ree ' collision integrals tore than double 

nrieger s values were obtained at low T*. At hi ah T* t>,a j.» 0 

lessens because the effect of polarity Lreases^A study of 

limits 6 *? 10 T -l ab1 "' ” lethod used b y Krieger indicated an error in the 
limits of an integration, and that the transformation to goniometric 
variables is unfeasible when 6 is greater than zero. D^ls 2e 
gxven in appendix B. The only other work for 5 greater than zSo is 

™iues 1 S ed s eal0Ul8 “° nS by Uh include* 1 negative 

Se ^sSts^f^S 1 ^ 5031 * 1 ’' 6 - Bl6lr r6SUltS ^ «th 


DETERMINATION OF PARAMETERS 

In order to determine the constants a e/k and b f*nr vn-rir.no 
molecuies Krieger assumed the molecules interred as poirt dSSJs of 
maximum attraction. He then selected two experimental vUcosufSta 

points for each molecule, and used equations (6) and (7) and his o(2,2)* 

a 6 /g ^ nectlon with the experimental data to determine the constants. 

HlS ** ’ table extends over a range of T* fmm i tn ri? „ n 

range of 5 from 0 to 2 at intervals of 0. 25. Of 12 molecules tested 

wa er had the highest 6, with a value of 2.33. In general the constants 
seem reasonable. ° di ^ne constants 

^ Procedure to find the constants using the present fl( 2 ; 2 )* 
table was the same as Krieger’s method, except that a least squares fit 
ec e experimental viscosity data was used to determine the best 
overall constants for a molecule. Table III glves the c™s?LS (a* 

e/k, 6) obtained using the present n( 2 ; 2 )* table and experimental 
dipole moments (p). Table IV gives the experimental viscosities used 
to Obtain the constants of table HI. It also contains Sscosities 

calculated with these constants and 2 )* values of this naner ms. 

agreement between experimental and calculated data is goS Se ^erS 



12 


and 


deviation for all molecules is 0.5 percent, the largest average deviation 
being 1. 2 percent for H 2 S. The explanation for the relatively large 
deviation for HgS is that the experimental data are not smooth. The 
constants, a and e/k, are different from those of nonpolar molecules; 
a is larger and e/k is smaller. This becomes more pronounced with 

increasing 5 values. 

No satisfactory results were obtained for more highly polar molecules 
such as NH* or HoO using this least-squares technique. Independent hand 
calculations verified the computer results and indicated that extending 
the SSes to larger values of 5 would not help. Since the contribu- 
tion of the dipole -dipole interaction term to the Stockmayer potential 
is small for slightly polar molecules, and becomes ^Portant for highly 
polar molecules, it appears that assuming g.0 1 ,e 2 ,cp) equals +2 at al 
times is inadequate. 

However, another possible method for obtaining the constants a 
6 would be to treat & as a third parameter, independent of 
and e. This would mean that no specific reLative orientation is as- 
sumed during the encounter. Then, knowing tie dxvole moment of th 
molecule, it would be possible to calculate an effective g{6 1 ,e 2 ,<p) to 
each molecule by equation (4). 

Hornig (ref. ll) suggests using a combination of ^ree types of 
interactions. Two are for resonant collisions, where the fir 
g(e 1 ,e 2 ,cp) equal to some positive number between 0 and 2, and the second 

is -g(0 1 ,0 2 ,cp). The magnitude of g( Sq, e 2 ,<P) is calculated from a. 

knowledge of the internal quantum numbers of the molecules. The third 
TZ of interaction is a nonresonant collision where the r 5 term 
disappears. Mathematically, this latter interaction is identical to 
the Lennard-Jones equation for the interaction of nonpolar molecules. 
Therefore, according to Hornig, the effect of the r”3 term on the 
potential can be attractive, repulsive, or nonexistent depending upon 
the type of interaction between the two molecules. 

Using Hornig' s approach, an effective g{9i,92,ty) for 811 interaction 

could then be calculated as a weighted average of the three types of 
interactions, where the weighting factor would depend upon the frequency 
of each type of collision. When collision integrals for negativ 5 
values become available, it will be possible to try this approach. 

In summary, it is concluded that the ebility to estimate transport 
properties of polar gases is still in doubt. However, two P° s ^ e 
approaches to the solution of the problem have been mentioned that may 
eventually resolve the situation. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland^ Ohio, August 12; I960 
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APPENDIX A 
SYMBOLS 


A l/a* 

a distance of closest approach of colliding molecules 

aQ^a-^ constants in empirical equation (35) for curve-fitting 

6^ against b*2 

a* a/ a 


b distance of closest approach of colliding molecules in 

absence of any interactions 


b* 

b/cr 


*4 

(l 

lover limit for analytic integration of Q v ' 
in vicinity of orbiting 

integrals 

b N 

upper limit for analytic integration of 
in vicinity of orbiting 

integrals 


numerical integration limit for integrals 

b 0 

value of b* for which orbiting occurs 


Ci(x) 

'/’ as ^ dt 


b l 

first approximation to coefficient of diffusion 


E(r) 

interaction potential 


E* 

E(r*)/e 


g 

relative velocity between molecules at infinite 
before colliding 

separation 

g* 

(mg 2 / Ic) 1 / 2 



numerical integration limit for ft(U s )* integrals 
value of g* such that orbiting does not occur for all 

S*>gQ 
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k 

M 

m 

P 

4 (l) (g,S) 

r 

* 

r 

Si(x) 

T 

rp* 

r 

6 

€ 


^1 

e 

e m 

Q±jQ2’V 


e 

p 

p 

a 

x 


Boltzmann ’ s constant 
molecular weight 
molecular mass 
pressure 

collision cross section 




rrigid sphere 

intermole cular separation Between colliding molecules 


r / o 

/ 

0 


sin t 
t 


dt 


temperature 

kT/e 

(mg 2 /4kT) 1 / 2 

parameter in modified Stockmayer potential 


maximum energy of attraction between colliding molecules 
in absence of dipole forces 


first approximation to coefficient of viscosity 

1 (* - X) 

2 v 

angle 9 at infinite separation of colliding molecules 

angles describing relative orientation of two point 
dipoles 

first derivative of 9 with respect to time 
dipole moment 

l/r* 

collision diameter 

angle of deflection in bimolecular collision 


t?d 
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fl(^ s )(T,5) 
a(i, s ) *(t*, s) 


collision integral 

n(^s)/ fl (z_,s) 

' rigid, sphere 
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APPENDIX B 

INVALIDITY IN THE TRANSFORMATION TO GONIOMETRIC VARIABLES 

In reference 3, the following statement is made on page 18: "The 

value of t>* = 0 (central collision), for which 2/a* 6 = 1 + (l + g* 2 ) 1 / 2 , 
corresponds to the value P = 0. " It can be shown that this is a true 
statement if, and only if, 5=0. To accomplish this, equations (26), 
(29), and (30) from reference 3, corresponding here to equations (Bl), 
(B2)^ and (B3), respectively, are used: 


+ - a * -6 - 6a *' 3 ) ■ 1 

(Bl) 

g* = cot r (o < r < 

(B2) 

2a " 6 - 1 + iS"r <° £ % + r) 

(B3) 

where equations (B2) and (B3) are the defining equations 
formation to the goniometric variables (3 and T- 

for the trans- 

Xf X)* = 0 (the lowef limit of the cross- section integrals equa- 

tion (Bl) becomes 


a*- 12 _ a* -6 - 5a*" 3 = 

(B4) 

Let a.Q be the one positive real root of^this equation, 
equation (B3) the 6 corresponding to b* = 0 is 

f“> \ “"I 

Then from 

P S arc cos N 2a*" 6 - ljsia TJ 


or 

P = arc cosU* ^ - l)/(g* 2 + l)^j 

(B5) 

since sin T = (g* 2 + l)" 1 / 2 from equation (12). 


If 6 ss 0^, equation (B4) becomes 


a*" 12 - a*" 6 = 

(B6) 


-791 
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The solution of this equation leads to the result 

a *-6 _ 1 i (l + g*2)l/2 
2 

and^ since aQ is positive, 

2a£~6 = 1 + (l + g *2)l/2 

If this quantity is substituted in equation (B5), the result is 

3 = arc cos 1 


or 


(B7) 


3 = 0 


When 8=0, 3 = 0 corresponding to b* = 0. 

Next, assume 3=0. Then from equation (B5), 

(2a£-6 - l)/(g*2 + !)l/2 = x 

Solving for a*-6 and substituting in equation (B4) give 



which, after combining like terms, results 


m 




;* 2 + l ) 1 / 2 + i j 1//2 


= 0 


(B8) 


(B9) 


This implies either 8=0 or (g*2 + i)l/2 + i = 0 . Slnce the latter 

® y ^ n6Ver Z6r ° for all / eal «*, & - 0- Thus it has been shown 

that 3=0 corresponds to b* = 0 if, and only if, 8=0. In ref- 
erence 3 3 - 0 is used for the lower limit of the cross-section inte- 

gral for all 8. This is therefore incorrect. 


However, even if the correct lower limit (eq. (B5)) for 3 had 
been used the transformation to 3 (eq. (B3)) does not always lead to 
a real value for 3 corresponding to b* = 0. An example will suffice 
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to illustrate this fact. Take 6=1 and g* = V 4 ® 921,4 calculate (3. 

Substituting in these values for g* and 6 and letting a* 3 = C, 
equation (B4) becomes 

C 4 - C 2 - C - 10 = 0 (BIO) 

The solution of equation (BIO) for the one real, positive root is 
C = 2. Therefore, 2C 2 = 2a*" 6 = 8. Substituting this value in equa- 
tion (B5 ) leads to the result 

„ 7 

P = acre cos — — : 

V 4 i 

Thus in this particular example, no real P corresponds to b =0. 
Therefore, the transformation to the goniometric variable p by equa- 
tion (B3 ) is not valid for 5^0. 
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Figure 5. - Angle in vicinity of orbiting. 




